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Abstract 

The focus of this paper is on the analysis of a semi-infinite crack lying along a perfect interface 
in a piezoelectric bimaterial with arbitrary loading on the crack faces. Making use of the extended 
Stroh formalism for piezoelectric materials combined with Riemann-Hilbert formulation, general ex¬ 
pressions are obtained for both symmetric and skew-symmetric weight functions associate with plane 
crack problems at the interface between dissimilar anisotropic piezoelectric media. The effect of the 
coupled electrical fields is incorporated in the derived original expressions for the weight function ma¬ 
trices. These matrices are used together with Betti’s reciprocity identity in order to obtain singular 
integral equations relating the extended displacement and traction fields to the loading acting on the 
crack faces. In order to study the variation of the piezoelectric effect, two different poling directions 
are considered. Examples are shown for both poling directions with a number of mechanical and 
electrical loadings applied to the crack faces. 

Keywords: Interfacial crack. Piezoelectric bimaterial, Weight function, Betti’s identity. Singular 
integral equations. 


1 Introduction 


Fracture in piezoelectric materials is an area of great interest due to their many applications in industry, 
for example, in piezoelectric insulators and actuators [H, 0|. The problem of a static semi-infinite interfa¬ 
cial crack between dissimilar anisotropic piezoelectric materials under symmetric loading conditions has 
been studied in Suo et al. Q using an approach based on the Stroh formalism Q and Riemann-Hilbert 
formulation. As an alternative to this method, singular integral formulations for two-dimensional inter¬ 
facial crack problems in mezoelectric bimaterials have been derived by means of approaches based on 
Green’s function method 0. Applying this procedure, the displacements and the stresses are defined by 
integral relations involving the Green’s functions, for which explicit expressions are required. Although 
Green’s functions for several crack problems in piezoelectric bimaterials have been derived , their 
utilisation in evaluating physical displacements and stress fields on the crack faces requires challenging 
numerical estimation of integrals for which convergence should be asserted carefully. Moreover, both the 
complex variable formulation proposed by Suo et al. and the approaches based on Green’s function 
method work when the tractions applied on the discontinuity surface are symmetric, but not in the case 
of asymmetric loading acting on the crack faces. 

In this paper, we illustrate a general procedure for studying plane interfacial crack problems in 
anisotropic piezoelectric bimaterials in the presence of a general non-symmetric mechanical and elec¬ 
trostatic loading distribution acting on the crack faces. The approach developed in Piccolroaz et al. 
0 and Piccolroaz and Mishuris 0, based on weight function theory and Betti’s reciprocity theorem, 
is extended to the case of an interfacial crack in a piezoelectric bimaterial containing a perfect inter¬ 
face. Applying this method, the crack problem is also formulated in terms of singular integral equations 
avoiding the use of Green’s functions and the resulting challenging computations. 

In fracture mechanics the notion of weight function, defined as the stress intensity factor associated 
to a point load acting on the crack faces, was originally introduced by Bueckner [m . fl^ and Rice 
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M- This weight functions concept have been extended for studying crack problems in piezoelectric 
materials by McMeeking and Ricoeur 141 and Ma and Chen [l^ . An alternative formulation where the 
weight functions are defined as the singular displacement field of the homogeneous traction free problem 
was proposed by Willis and Movchan [^. Recently, this weight functions definition has been used in 
the derivation of stress intensity factors for both static and dynamic crack problems in isotropic and 
anisotropic bimaterials iEQ and in thermodiffusive elastic media [l^ . These weight functions 
have also been used in the derivation of singular integral equations relating interfacial tractions and crack 
displacement to the applied loadings on the crack faces [9|,l20|,[2l| . In this paper, this method is generalised 


in order to study interfacial cracks in piezoelectric bimaterials. Explicit weight functions are derived for 
the class of transversely isotropic piezoelectric bimaterials considering two different poling directions 22| . 
Then these weight function matrices are used together with the generalised Betti’s identity to derive 
singular integral equations relating the tractions and displacements to the electric fields introduced by 
the piezoelectric effect. Finally, some simple illustrative examples of the application of the obtained 
singular integral identities are reported. 

Section 2 of the paper introduces the problem geometry and a number of preliminary results required 
for further analysis. In Section 3 we derive weight functions for transversely isotropic piezoelectric bima¬ 
terials considering two different poling directions using the definition introduced by Willis and Movchan 
0. The problem is reduced to only contain the fields affected by the piezoelectric effect. In Section 4 
the obtained weight function matrices are used together with Betti’s identity in order to derive singular 
integral equations relating the physical fields to the applied mechanical and electrical loadings on the 
crack faces. In Section 5 we present a number of examples for a variety of mechanical and electrical load¬ 
ings. For one of the considered examples we also present the results from finite element computations 
and compare their accuracy to those obtained from our singular integral equations. 


2 Problem formulation and preliminary results 


In this section we introduce the mathematical model used for the remainder of the paper. Some pre¬ 
liminary results concerning two-dimensional interfacial cracks in anisotropic piezoelectric bimaterials are 
reported. 

We consider a semi-infinite crack lying along a perfect interface between two dissimilar piezoelectric 
half-planes, referred to as materials I and II. The crack occupies the region {xi < 0, X 2 = 0}, as illustrated 
in Figure [T] The perfect interface conditions in a piezoelectric bimaterial are continuity of displacements, 
traction, electric potential and the electric charge. The loading along the crack faces, for xi < 0, is known 
and given by the functions 

p^(xi) = cr2j(xi,0=^), for j = 1,2,3, p^{xi) = (I) 


where aij and Di represent tractions and electrical displacements respectively. 

Expressions for the stress fields and displacements for a plane semi-infinite interfacial crack between 
dissimilar piezoelectric media have been derived by Suo et al. Q, using an approach based on Stroh 
formalism and Riemann-Hilbert formulation. This method, which is applied in the paper in order to 
derive explicit weight functions for piezoelectric bimaterials, is summarised in Section l2.Il 

In Section 12.21 the definition of weight functions proposed in Piccolroaz et al. 23| and Morini et al. 


I7| is extended to the case of interfacial cracks between dissimilar piezoelectric materials. Finally, in 


Section nn the expression for the Betti’s formula generalized to piezoelectric media by Hadjesfandiari 
is reported. Further in the text, weight functions and Betti’s formula are used for formulating the 
crack problem shown in Fig. [T]in terms of singular integral equations. 


2.1 Stroh formalism for piezoelectric materials 


Under the conditions of static deformations, the governing equations for a linear and generally anisotropic 
piezoelectric body are [2^, [26| : 

D^,i = 0 , ( 2 ) 


^ij,i — bi 


where body forces are assumed to be zero. The strain e and electric field E, are defined by the gradients 




E. = 


(3) 


2 










C/,w/,e/ 


vt 


Pt 


Pt 


X2 


Et 




crack 


interface 

->- 

Xi 


P2 


E 2 




Cii,uJii,ejj 


Figure 1: A semi-infinite crack along an interface between two dissimilar piezoelectric materials subject 
to the state of generalised plane strain and short circuit (ef = Ef = 0 ) 


where Ui are the components of the displacement field and (/> is the electric potential. 

The constitutive relations for anisotropic piezoelectric materials are given by Wang and Zhou [ 2 ^ 

(Tij — CijfsSrs (^sji ^st Di — CJis^ s H” ^irs^rs j (4) 

where the tensors C, lo and e represent the stiffness, permittivity and piezoelectricity respectively. Com¬ 
bining equations ([U, ([3]) and (|1]) the following relationships are obtained 


{CijrsUr -f ,si — 0, ( ^is4^ “b Girs'^r') ,si — 0- 


(5) 


Considering the semi-infinite interfacial crack problem shown in Fig. [1] and assuming for both the 
upper and lower piezoelectric half-planes the state of generalised plane strain and short circuit (£3 = 0 
and A 3 = 0), a solution for equations (O is derived by applying the procedure described by Suo et al. 
Q. The solution is sought in the form {ur,4>)'^ — ^/(z) where z = xi + px 2 - This yields the following 
eigenvalue problem 

[Q-bp(R-l-R^)+p2T]A = 0 (6) 

For general, anisotropic piezoelectric media the matrices have the following form: 


/Cn 

Ci6 

Ci 5 

eii \ 


^C'le 

C 12 

Ci 4 

616 ^ 


^(^66 

^26 

C46 

626 \ 

Ci6 

Cee 

C56 

ei 6 

, R = 

Cqq 

C26 

C46 

612 

, T = 

C26 

C22 

C24 

622 

C'ls 

Cse 

C55 

ei 5 


Cse 

^25 

C45 

614 


C46 

C24 

C44 

624 

Veil 

ei 6 

ei 5 

-Wii/ 


Ve2i 

626 

625 

-W12/ 


V626 

622 

624 

— W22/ 


The eight eigenvalues of are found to be four pairs of complex conjugate 

For the remainder of this paper u shall be used to the extended vector containing both the physical 
displacement and electric potential given by u = (ui, U 2 , U 3 , (^)^. The extended traction vector, t, is also 
introduced, given by t = ((T 2 i, £* 2 )^. It was shown in Suo et al. Q that the extended displacements and 
traction are given by 

4 4 

u = 2Re^A^/^(z^), t = 2Re^L^/;(z^). (7) 

fl—1 fl—1 


where 


3 

^ (^lj 2 P/j.C 2 j 2 ; for j = 1,2,3, 

r—1 
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Li^ = ^^[(e2T-l + P/i,e2r2)^r/i] — (wi2 + P/i,W22)^4/i ■ 
r—1 

The values of are given by = xi + Pfj,X 2 , where the eigenvalues, p^, are taken to be those with 
positive imaginary part and the matrix are the corresponding eigenvectors. The surface admittance 
tensor B = iAL~^ is introduced and for further analysis the following bimaterial matrices are defined 

H = B 7 + B'//, W = B/-B//. ( 8 ) 

Assuming vanishing tractions and electrical displacement along the crack faces for xi < 0, and using 
expressions in Suo et al. Q the following Riemann-Hilbert problem was found to exist along the 
negative xi- axis 

h^(a;i) + H“^Hh“(a;i) = 0, —oo < a;i < 0. (9) 

A solution is found in the form h( 2 ;) = where the branch cut is situated along the negative real 

axis (which is the crack line). Inserting this solution into equation ([2]) yields the eigenvalue problem 

Hw = (10) 

Four pairs of eigenvectors and eigenvalues are found from (US. They are 

(e,w), (-e,w), (-ik,W 3 ), (zk,W 4 ). ( 11 ) 

Imposing the continuity of tractions and electrical displacement along the interface ahead of the crack 
tip, the following formula for determining the extended traction vector t is derived [^: 

h^(a;i) + H^^Hh“(a;i) = t(a;i), 0 < xi < oo. ( 12 ) 

The extended traction vector, derived by means of equation (ini, is therefore given by 

t(a;i) = {2ttxi)~^ [Kx'^i'w + + K^x^ws + K^x'^'^'W/^, (13) 

where K = Ki + 1 X 2 - The stress intensity factor vector is given by K = {K, K, A 3 , . 

The extended displacement jump across the crack, given by |u]|(a;i) = u+(xi) — u“(a;i), was found 
to be 


H(a^i) = (H + H) 


(—cci) A(—xi)*'^w K{—xi) 

27r [(1 + 2ie) coshTre (1 — 2ie) cosh-Tre 

K:i{-xi)''^^3 ^ A4(-xi)~"^W4 

(l + 2K)cos7rK (1 — 2 k) cos7rK_ 


+ 


(14) 


Using these expressions for the traction and displacement jump the following equation can be used to 
find the energy release rate at the crack tip 


t^(A-r)[u](r)dr, 

for an arbitrary value of v. Using expressions m 9(iici (jl4p it c&n be shown thnt 

Wg" (H + H)W4 


G = ^d2h^iAf 

4 cosh Tre 


4 COS^ TTK 


-A 3 A 4 . 


(15) 


(16) 


2.2 Weight functions for piezoelectric bimaterials: definition 

Weight functions, introduced in Bueckner [HI [I^, are functions whose weighted integrals can be used 
in the derivation of important fracture parameters, such as the stress intensity factors. For elastic 
materials Willis and Movchan [l^ introduced a weight function given by the singular displacement field 
corresponding to a homogeneous, traction-free problem similar to Fig. [T] with the crack occupying the 
region a:i > 0 and the perfect interface lying along the region Xi < 0. This concept of a weight function has 
been adopted in the works of Piccolroaz et al. and Morini et al. 0 for studying interfacial cracks 
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in isotropic and anisotropic bimaterials, respectively. Here, this approach is extended to piezoelectric 
materials. 

The weight function for piezoelectric materials is given by the extended singular displacement field U 
incorporating both displacement and electric potential. The symmetric and skew-symmetric parts of the 
weight function across the plane X 2 = 0 are given by 

[Ul(xi) = U(a:i,0+)-U(xi,0-), (17) 

(U)(xi) = i [U(a;i,0+) +U(ii,0-)] . (18) 

To satisfy the perfect interface conditions it is clear that |U] = 0 for xi < 0. 

The extended traction field corresponding to the extended displacement U, is denoted S. The fol¬ 
lowing Riemann-Hilbert problem is found along the positive portion of the xi-axis 

h+(a;i)-f H-iHh-(a;) = 0, 0 < xi < oo. (19) 

A solution for h{z) is now sought in the form h = The branch cut of h is situated along the 

positive part of the xiaxis. Inserting this solution into (IT^ yields the following eigenvalue problem 


Hv = (20) 

It is immediately clear by comparing equation m to (pH)) that V = w, V 3 = W 4 and V 4 = W 3 . 

Along the negative part of the real axis S is given by 

h+(xi)-f H-iHh-(x) = S(xi), -oo<xi< 0 . ( 21 ) 

Therefore the extended traction vector corresponding to the weight function U is given by 

(-X,)-3/2 

S(a:i) = — ^ [C(-xi)^‘^w-b C(-xi) + Csi-xi) '‘W 3 -f C' 4 (-a;i)”w 4 ] , (22) 

2 v Stt 

where C = Ci + iC 2 , C 3 and C 4 are constants defined in the same manner as the stress intensity factors 
for the physical problem. 

Let us introduce the Fourier transform of a generic function / with respect to the variable xi , defined 
as follows: 

/ +00 1 ^- 1-00 

/(xi)e*«"M:ri, f{x,)=J^-^[f{0] = ^ (23) 

-00 J —00 

For anisotropic materials it was shown in Morini et al. 0 that the Fourier transforms of the symmetric 
and skew-symmetric parts of the weight functions are related to S in the following manner 


IU ]+(6 = ^(*sign(e)Im(H) - Re(H))E-(0, (24) 

(U)(C) = ^(*sign(C)Im(W) - Re(W))E-(0, (25) 

As the method used in Morini et al. 0| was for general matrices H and W, it is immediately clear 
that these results also hold for piezoelectric materials. 


2.3 The generalised Betti formula 


In this part of the paper we consider the Betti identity in the context of a semi-infinite crack in a 
piezoelectric bimaterial. Originally used to relate two different sets of displacement and traction fields 
satisfying the equilibrium equations [l^, [ 2 ^ , this approach was extended to study piezoelectric materials 
(with electric potential and electric displacement) by Hadjesfandiari 2^. The derivation of Betti formula 
for piezoelectric solids is briefly reported in Appendix This integral formula is now used to form a 
relationship between the physical fields and the weight function introduced in the previous part of the 
paper. 
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Applying the generalised Betti formula derived in Hadjesfandiari [2J| to the semi-circular domain 
occupying the upper half plane in the model illustrated in Fig. [U the following integral relation is found 


j [RU(a :2 — xi, O'*') • t(a:i, O'*') — RS(xi — xi, 0+) • u(xi, 0'*')] dxi = 0, 

J 2:2=0+ 


(26) 


where u = {ui,U 2 ,U 3 ,(j))'^ and t = (cr 2 i, cr 22 , 023 , £> 2 )^ (see the Appendix for details) and R is given by 


/-I 0 0 0 \ 

0 10 0 
0 0-10 
V 0 0 0 -1/ 


The equivalent equation for a semi-circular domain in the lower half-plane gives 

/ [RU(x^ — xi, 0“) • t(xi,0~) — RS(x^ — xi,0~) • u(xi, 0“)] dxi = 0. (27) 

J X2—0~ 

Subtracting equation (1271) from (1^51) yields the following relationship 

RIU] * t(+) - RS(-) * H = -RIU] * (p) - R(U) * [pl, (28) 

where * represents the convolution with respect to xi and is used to represent the restriction of a 
function to the positive or negative portion of the xi-axis respectively. It can be easily deduced that 
in equation (1281) the contribution to the generalised traction vector defined on the negative semi-axis 
xi < 0 is given by the loading functions = (o27 \022 \ 023 \ = {Pi,P 2 ,P 3 ,P 4 )'^ = P, and the 

symmetrical and skew-symmetrical part of the load, respectively (p) and |p], are defined as follows [^: 

(P) = ^(p-^+p-), Ip1=P+-P~- (29) 

Applying the Fourier transform to (pS)) gives the following relationship 

lUf Rt+ - (S^)^RIuI = -lUf R(p) - (U)^R[pl. (30) 

In the next Sections, explicit expressions for the weight function matrices (1241) and dSH]) are derived and 
used together with the the generalised Betti identity dim for formulating the considered interface crack 
problem in terms of singular integral equations. Since the bimaterial matrices H and W involved in the 
weight functions (HI and (1251) depend on the surface admittance tensors of both piezoelectric half-planes 
(see definition (©), in order to derive explicit expressions for these matrices the solution of the Stroh’s 
eigenvalue problem ([5]) is needed. In the general fully anisotropic case, this eigenvalue problem must be 
solved numerically. Nevertheless, exact algebraic expressions of Stroh’s eigenvalues and eigenvectors have 
been obtained for the class of transversely isotropic piezoelectric materials in [i,[2i,[2i| and [1^. This 
class of materials has practical significance, because many poled ceramics that are actually in use fall into 
this category. The Stroh matrices and surface admittance tensors for transversely isotropic piezoelectric 
materials with poling direction parallel to X 2 and X 3 axes are reported in the Appendix. Further in the 
text, these results will be used together with expressions (El and (1251) for deriving explicit weight function 
matrices corresponding to interfacial cracks between dissimilar transversely isotropic piezoceramics. 


3 Weight functions 

Piezoelectric materials occupying both lower and upper half-planes in Figure [I] are assumed to be trans¬ 
versely isotropic. For simplicity, poling direction is assumed to be parallel to the X 2 and X 3 axes, respec¬ 
tively. Using eigenvalue matrices and surface admittance tensors in the forms reported in the Appendix, 
explicit weight functions are derived for both these cases. 
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3.1 Poling direction parallel to the 2 : 2 —axis 


Poling direction directed along the X 2 —axis is assumed for both upper and lower piezoelectric half-planes. 
Considering the geometry of the model shown in Figure[T] it is easy to observe that in this case the poling 
direction is perpendicular to the crack plane. Under these conditions the Mode III component of the 
solution decouple from Modes I and II and the piezoelectric effect 22, MM , which further in the text we 
will refer to as Mode IV. This means that the antiplane tractions and displacement have no dependency 
on the electric field and therefore behave in the same way as they would in an elastic material with no 
piezoelectric effect. The stiffness, permittivity and piezoelectric tensors corresponding to this case are 
reported in the Appendix together with explicit forms of the matrices involved in the decoupled part of 
the eigenvalue problem ©■ 

The remainder of this section of the paper considers only the in-plane components and electrical effects. 
That is: u = (iti, U 2 , 4’)'^ t = ((T21, o'22, £^2)^- The surface admittance tensor B then becomes 



/ 

Bii 

iBi 2 

iBi 4 

B = 


-iBi 2 

B22 

B24 


V 

-iBi 4 

B24 

B44 


(31) 


The expressions for the components of the matrix B, found by Hwu for the two-dimensional state of 
generalized plane strain and open circuit, are quoted in the Appendix of the paper. With an expression 
for B it is now possible to construct the bimaterial matrices required. The bimaterial matrices H and 
W can be written as 



( H41 

iHi2 

lHi 4 \ 


( VFii 

iWi2 

iWi 4 \ 

H = 

-iHi2 

H22 

H24 

, W = 

-iWi 2 

IU22 

W24 


\-iHi 4 

H24 

H44) 


\-iWi 4 

IU24 

W44) 


Note that explicit expressions for the components of matrices are also reported in the Appendix. 

Knowing the structure of the bimaterial matrix H it is now possible to find expressions for the traction 
field, E using the eigenvalue problem (HOD. From ([321(1) it is only necessary to find three sets of eigenvalues 
and eigenvectors. They have the form 


(e,w), (-e,w), (iK,W 4 ). 

As the Mode III components of the solutions have decoupled and behave purely elastically it is expected 
that At = 0. Therefore it is expected that the eigenvalues are given by two non-zero real valued numbers, 
with the same magnitude but differing in sign, and 0. With these particular eigenvalues and eigenvectors 
the expression for E is given by 

S(a:i) =- -= — [C{-xiy''v/+ C{-xi) -|-C' 4 W 4 ] . (33) 

2\/27r 

To find the eigenvalues from m the following equation must be solved 

IIH-e^’^'^HII = 0. (34) 

Substituting ([22) (1) (IMl) the following equation is derived 

(1 - - e^^y^Hn(H22H44 - HI 4 } - (1 + e‘^^y\Hl^H44 + HI 4 H 24 - 2 H 42 H 44 H 24 )] = 0. (35) 

As expected solving the equation 1 — = 0 yields the eigenvalue 0. The other eigenvalues are given by 

where 

= 41 = H44{H22H44 - HI 4 ). B = 2{HI^H44 + HI^H24 " 2H42Hi4H24). 
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Using these eigenvalues it is possible to find expressions for the eigenvectors w and W 4 . The expressions 
chosen here are made for notational convenience. The expression for W 4 is 


1 


/ 0 


W 4 = 


£fl4 

V-^i 2 y 

There are three possible expressions for the eigenvector w. They are 


(37) 


1 

w — - 
2 


f-iPiHh - 

H 44 H 12 — H 14 H 24 
\ H 14 H 22 — H 24 H 12 


1 

2 


^—il3{Hi4H22 — H 42 H 24 ) 
l3^HiiH24 — H 12 H 14 
11 H 22 


Hu H 2 


1 

2 


/-Z/3(i7l4i724 - H42H44f 
P^HuH44 - H^4 
\ H 12 H 14 — f3‘^HiiH24 


(38) 

For the remainder of this paper the first representation of w from equation (1381) shall be used. 

Using (1551) it is possible, using the method described in Piccolroaz et al. 0, to construct three 
independent traction vectors using the following three cases: 

1 . Cl = 1,C2 = C4 = 0 , 

2 . (^2 = 1, Cl =6-4 = 0, 

3. C 4 = l,Ci = C 2 = 0 . 

Using dSHD the three traction vectors obtained are 






. 3/2 / imh - H22H44)[i-xiy^ - (-xi)--] 

^ {H44H12 - Hi 4 H 24 )[{-Xiy^ + {-xi)-^y 

\{Hi 4H22 - H24Hi2)[{-Xiy^ + (-Xi)—] 


2 v^ 

(-xi)-^/^ 


/ -/3(i7|4 - H22H44)[i-Xiy^ + (-Xi)--] 

*(i744-ffl2 - Hi4H24)[{-Xiy^ - {-Xi)-^y 
V(i7l4-ff22 - H24Hi2)[i-Xiy^ - (-Xi)-*^] 

/ 0 

Hi4 

\-Hi2j 


S^xi) = 


-3/2 


(-^ 1 ) 

2 v^ 


(39) 


(40) 


(41) 


Here, a superscript 4 has been used instead of 3 in equation (HD) so as not to confuse this with the Mode 
III components which have already been decoupled. 

In order to calculate explicit expressions for |U]|^ and (U) it is necessary to find the Fourier transforms 
of dMl), (gOl) and (HU) 




(1 + 4e2)0r 


= 


V2ii^y^^ 
(1 + 


^ iy{H24 - H22H44) [(-i - ze)F(i + ie){iy) - (-5 + Je)r(i - ie){i^yy 

{H44H12 - H14H24) [(-^ - *e)r(i + + (-5 + *e)r(i - ze)(iO"] 

\{Hi 4 H 22 - H24H12) [(-i - ie)r(i + ze)(fO— + (-i + ^e)F(i - ze)(zO*l 

^ -P{H24 - H22H44) [(-^ - ie)r(i + + {-\ + - ie)(i^)*'-] 

i{H 44 Hi 2 - H14H24) [(-i - ie)r(i + ie)(zC)-*^ - (-i + ie)T{\ - ie)(iO"] 
\i{Hi 4 H 22 - H24H12) [(-i - ie)r(i +ie)(zC)-*^ - (-i +ie)r(i - ze)(ioH, 


(42) 


= 


v /2 


( ^ 

-Hi4 

V Hi2 


With these expressions it is now possible to use a 3x3 matrix whose columns are the three linearly 
independent traction vectors found along with equations (1551) and (1551) to find expressions for JU] and 


















3.2 Poling direction parallel to the X 3 —axis 

Observing Figure [TJ it can be noted that in the case where both the upper and lower piezoelectric half¬ 
planes are assumed to be poled along the xa—axis, the poling axis coincides with the crack front. For 
this particular case it is possible to decouple the Mode I and Mode II components of the displacement 
and stress fields from the Mode III fields and piezoelectric effects on the material 2^ 221. This means 
that the in-plane fields will behave similarly to those for purely elastic materials with no piezoelectric 
behaviour. Also for this case, the explicit form for the stiffness, permittivity and piezoelectric tensors are 
reported in the Appendix. 

In the remainder of this section only the out-of-plane and piezoelectric components are considered. 
That is: u = {u 3 ,(j))'^ and t = {(J 23 , D 2 )"’"■ The surface admittance tensor then becomes 


B = 


B33 B34 


B34 B441 

Consequently, the bimaterial matrices H and W can be computed and have the form 


H = 


H 33 H34\ _ f ^3^33 JH 34 \ 

^H34 H44)' ~\iH34 S4H44)' 


(45) 


(46) 


where: 


Hap = [Bap\i + [Bq,/?]//, for a, /3 = 3,4, 


5a = 


[Bckq]/ [Bcto-]// 

[Baa\l + [Baa]ll ’ 


for a = 3,4, 


7 = 


[^34]/ — [B34]ii 

H34 


Explicit expressions for the components of Bap, are given in the Appendix of the paper. 

In order to obtain the weight functions for the materials considered here the Riemann-Hilbert problem 
([5]) must again be considered. For this case the bimaterial matrix H has no imaginary part, and then 
substituting expression (021) (1) into equation dH]) we get 


h^(a;i)-|-h (x) = 0, —00 < xi < 0. 


(47) 


For this special case it was shown in Suo et al. Q that the extended traction along the interface and 
displacement jump across the crack are given by 


t(a:i) = (27ra:i) = 


K3^ 

K4j 


for Xi > 0, 


(48) 


H(a;i) = 


2 (-xi) 


H 




for xi < 0. 


(49) 


Knowing the traction fields makes it possible to evaluate the weight function U and its corresponding 
traction S for these particular materials. The expression for S is 


S(xi) = 


(-xi) 



for Xi < 0. 


(50) 


The Fourier transforms of th symmetric and skew-symmetric parts of the weight function, |U] and (U), 
are once again given by equations ([M]) and (03. However, due to H, and therefore W being purely real 
the expressions simplify to 


[Ul+(e) = -^HS-(0, 


(51) 

(52) 
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where S is the 2x2 matrix consisting of two independent tractions. The linearly independent tractions 
are given by the case Cz = 1,0^ = Q and Ca = 0, (74 = 1 in equation dSO]). The results obtained are; 


S"(xi) = 








(53) 


The Fourier transforms of these vectors, which are clearly very important in deriving the weight functions, 
are given by 




, 0 I 




0 1 


(54) 


4 Integral identities 

In this Section the obtained weight function matrices are used together with the Betti identity (1501) 
to formulate the considered crack problem in terms of singular integral equations. Integral identities 
relating the applied loading to the resulting crack opening and traction ahead of the tip are derived for 
transversely isotropic piezoelectric materials in both the cases where poling direction is parallel to the X 2 
and xz axes. 


4.1 Poling direction parallel to the a; 2 —axis 


Considering the case where both upper and lower transversely isotropic piezoelectric half-spaces possess 
poling direction parallel to the a; 2 —axis (perpendicular to the crack plane), the in-plane fields and piezo¬ 
electric effect decouple from the antiplane displacement and traction. Consequently, the Betti formula 
still has the form 


[Uf Rt+ - S^R[ur = -[Uf R(p) - (U)^R[pl, 


(55) 


where |U] and (U) are given by expressions (l24ll and (1251) together with matrices (l32ll , and the rotational 
matrix R becomes 


R = 


/-I 

0 


VO 



Multiplying both sides of (1551) by R ^ |U]| ^ yields 


the following equation 


t+-N[ur =-(p)-M[p], 


(56) 


where M and N are given by 

M = R-i|Ur^(U)^R, N = R~i|U]|-^S^R. 


(57) 


Using (l24l) and (1251) full expressions for M and N can be found: 

M=^(M'-f*sign(C)M"), (58) 

N = M (N' + zsign(0N"), (59) 

where explicit expressions for D, M', M", N' and N" can be found in the Appendix of this paper. 

Taking the inverse Fourier transform of equation (I56L the following equations are found for xi < 0 
and xi > 0 respectively: 


•^."VoINIun = (p)(xi) + J--/<o[M[pl], XI < 0, (60) 

t(+)(a:i) -f J'-i>o[M|p]] = xi > 0. (61) 

The term involving t cancels for xi < 0 as it is only defined along the interface and the (p) does not 
appear for cci > 0 as it is only defined along the crack. In order to derive explicit expressions for equations 
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(EOl) and ((HTll . the inverse Fourier transform of isign(^)|u| , |^||u] and i'flu] are computed applying 
the convolutions theorem (see Appendix for details). 

The singular operator S and the orthogonal projectors T’± are defined 


Sip = -* ip{xi), 

TTXl 


(62) 


'P±ip 


ip{xi), ±xi > 0 , 
0 , otherwise. 


(63) 


Introducing the singular integral operator 5^®) = VSV- and the compact operator 5^'^^ = V+SV-, 
equations (15(71) and (1511) become 


A/'^^’®^ = (p)(xi) + M(^)[pl, x^<0, 

(64) 

\ XI >0. 

(65) 

The matrix operators and are given by 


7W(®) = ^ (M'-bM"5(®)) , 

( 66 ) 

Ar(®) = ^ (n'5(®) - N") , 

(67) 

2D 

( 68 ) 


(69) 


4.2 Poling direction parallel to the 0 : 3 —axis 

For the case where both upper and lower transversely isotropic piezoelectric half-spaces possess poling 
direction parallel to the ^3 axis, the weight functions consist of the 2 x 2 matrices m and dUl). The 
Betti identity (1301) then becomes a 2 x 2 matricial integral equation, where the rotational matrix R is 
given by 



Therefore, equation (13171) can be simplified further to give 



luf t+ - S^[ur = -[Uf (p) - (U)^[pl. 

(70) 

Multiplying both sides of equation (17(71) by JU] gives 



t+ - = -(p) - iur^(u)^[pr. 

(71) 

Using (I 5 TI) and (l52l) gives 

- z(oiur = -(p) - yipI”, 

(72) 

where 

Y=iH-iW, Z = -|e|H“i. 



Taking inverse Fourier transforms and methods seen in Piccolroaz and Mishuris @ and Morini et al. 
[13 gives the following singular integral equations 


-=-(p)(xi) - Y|p](a;i), forxi<0, (73) 

UXi 
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for Xi > 0 


t(xi) = -Q 


(c) 


dx\ 


(74) 


where and the integral operators 5^®) and 5^'^^ are defined in the same 

way as those introduced in previous Section. 

The integral identities dMl), (EH), (1231) and dUl) relate the mechanical and electrostatic loading applied 
on the crack faces to the corresponding crack opening and tractions ahead of the tip. The crack opening 
associated with an arbitrary mechanical or electrostatic loading can be derived by the inversion of the 
matricial operators Af*-®^ and Q^®^ in equations dSH) and 1721). Using the obtained crack opening functions 
in ESI and 17i)l. explicit expressions for the tractions ahead of the crack tip are yielded. Some simple 
illustrative examples of this procedure are reported in next Section. It is important to note that the 
solution of the obtained systems of singular integral equations El, El, Gl and 0, provides the 
crack opening displacements, the electric potential on the crack faces, the mechanical tractions and 
electric displacement ahead of the tip associate to general mechanical, electrostatic or electro-mechanical 
loading without any restriction concerning the geometry and the symmetry. In particular, the explicit 
evaluation of the skew-symmetric weight function matrices makes possible to consider the effects of skew- 
symmetric contributions to the loading. These effects cannot be accounted by means of other approaches 
available in literature [s, 3^, ^[l^, which are based on the assumption that the geometry of the applied 
loads is symmetric. 


5 Illustrative Examples 

In this Section we consider some examples of loadings for both poling directions and find solutions 
using the respective singular integral equations derived previously in the paper. Both mechanical and 
electrical configurations will be considered. Explicit expressions for crack opening and tractions ahead 
of the tip corresponding to both symmetrical and skew-symmetrical mechanical and electrostatic loading 
configurations are derived. The proposed illustrative cases show that the obtained integral identities 
represent a very useful tool for studying interfacial crack problems in piezoelectric bimaterials. To begin 
with we consider a symmetric distribution of point loadings when the poling direction is parallel to the 
X 2 -axis before considering both symmetric an asymmetric loading configurations for the piezoelectric 
bimaterial poled in the direction of the xa-axis. For the decoupled Mode III and IV example with 
symmetric loading we also present a comparison between the results from our singular integral equations 
and those obtained using finite element methods in COMSOL Multiphysics. 


5.1 Poling direction parallel to the X 2 -axis under symmetric mechanical load¬ 
ing 

We consider a symmetric system of two perpendicular point loads of varying magnitude on each crack 
faces acting in the opposite direction to their corresponding load on the opposite crack face at a distance 
a behind the crack. Mathematically these forces are represented as 


(P)(a:i) 


^—Fi5{xi + a)\ 
-F25{xi -fa) , 


V 0 


IpI = 0. 


(75) 


Under such a loading the singular integral equations used to find |u]| and t reduce to 

A/'(®)^^^l- = (p)(a:i), for xi < 0, (76) 

UX-y 

t(+)(xi) forxi>0. (77) 

UXy 

To simplify the problem we consider the set of bimaterials for which the matrix H from equation 
(l32U i^ has no imaginary part, that is i7i2 = Flu = 0. An example of when this may occur is when the 
upper and lower materials are the same. Under the matrix N" = 0 and therefore the integral equation 
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for xi < 0 becomes 


1 

D 


/Nil 0 0 

0 N22 -^24 I 

y 0 A^24 -^44 ) 


aiui( 


-) 


dxi 


^—Fib{xi + a)^ 
—F^bixi + a) 

V 0 


(78) 


From the system (1751) it is possible to obtain the following three equations for the derivatives of the 
displacements and electric potential 


^ = -FiDSixi + a), 


Sis) 




dxi 


dxi 


Sis) 


UXi 


= —F2D5{xi + a), 


dxi 


= 0 . 


(79) 

(80) 
(81) 


It is clear from these equations that the xi directed part of the solution for the example considered 
decouples from the component of the solution in the X 2 direction and the electrical effects. As a result of 
this we first proceed with solving for ui before proceeding to find expressions for U 2 and (j}. 

We begin by inverting the integral operator 5*-®^ in equation (1791) using the methods seen in Piccolroaz 
and Mishuris H and Morini et al. I 20 II 


dxi 


-) 


FiD 


Nun \ xi xi-r] 


Tj S{ri + a) , FiD j-a 
dr] = 


1 


Nun V xi xi + a 


(82) 


Using the expressions for D and Nn reported in the Appendix of this paper this expression simplifies to 

F,H 


11 


dxi 


a 1 
xi xi + a’ 


(83) 


which agrees with those results found in Piccolroaz and Mishuris Q and Morini et al. ( 2 ^ when a 
component of a displacement field decomles from all other components. Note here that the results differ 
from those in Piccolroaz and Mishuris [9j due to the anisotropy of the material and they differ from that 
in Morini et al. [ 2 ^ due to ui being decoupled in this paper whereas in that paper U 2 was separated from 
the rest of the solution. 

Integrating (l83l) gives the following expressions for the displacement jump along the crack 


|uil(a;i) = 


2Fii7i 


Xi 


■ arctanh w-, for — a < a;i < 0, 

a 


luiKxi) = 


2FiH 


11 


arctanh, / — 


Xl 


for Xl < —a. 


(84) 


(85) 


We now proceed to invert the operator in equations (1801) and m in the same manner. The resulting 
equations are therefore 


N 22 


dxi 


+ N 24 


dxi 


F 2 D I r 1 
n Y Xl Xl + a’’ 


N24 




■ T Af44 




= 0 . 


i9a;i dxi 

Solving these equations and simplifying gives the following expressions 


F2H22 r~d~ 1 

i 9 xi n Y Xl Xl + a’ 

dmi-i F2H24 1 

dxi n Y Xl Xl + a 


( 86 ) 

(87) 


( 88 ) 

(89) 
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Integrating these expressions gives the following expressions for the jump in displacement and potential 
over the crack 


) 


2F-2 


arctanh 


a \-H 24 J 


for — a < xi < 0, 


(90) 



2 F 2 


arctanh. / — 



for xi < —a. 


(91) 


Using equation ([77| the following expressions, for use in finding expressions for the interfacial traction 
and electric displacement, are obtained 




(Nn 

0 

0 \ 

022 

0 

N22 

fV24 

[02) 

1 0 

A^24 

N44J 


dxi 


The decoupled traction component can now be found: 


a 


(+) 

21 


(a^i) 


fVii 1 ^ El 

Dtt y_oo xi — rj drj tt \ xi xi + a 


(92) 


(93) 


Once again the obtained result is identical to that obtained for a decoupled field in anisotropic bimaterials 
[^ . with the only difference here arising from the difference in direction of the decoupled field. Using 
the same method the expressions for the coupled portion of the traction and electric displacement field 
are given as 

a22(a:i)(+) = i^^+^(xi) = 0. (94) 

It is seen that the mechanical part of the solution behaves identically to that in an anisotropic bimaterial 
and there is no electrical displacement component along the interface for any bimaterial with the given 
conditions. 


5.2 Poling direction parallel to the xs-axis 

5.2.1 Symmetric mechanical loading 

The loading considered here consists of a point load acting in opposite directions on each of the crack 
faces at a distance a from the crack tip. Mathematically this system of forces is represented using the 
Dirac delta distribution. The expressions for the symmetric and skew-symmetric parts of the extended 
loading are given by: 

(P)(a;i) = ^ : [pKa^i) = 0. (95) 

Inserting these expressions into equation (1231) gives the singular integral equation 


s(.») 


dxi 


= F 


H 33 

H 34 , 


S{xi + a). 


Inverting the operator 5*-®^ in a similar way to that seen previously in the paper gives 

_ F fH33\ llfS{r] + a) 


dx\ 


Integrating this equation gives the result 


TT \H34 J-oc \ xi Xi-r] 


-dry 


I 


H 33 


XiXi + a \ il 34 




2F 

= — arctanh ^ / — 


^ I H33 
a Vff34, 


for — a < xi < 0, 


(96) 


(97) 


(98) 
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I^al 


2F , , ^ 

= — arctann . - 



for xi < —a. 


(99) 


Making use of equation (1741) it is possible to obtain the expression for the extended traction vector, 
t, along the interface: 




H-1 

TT 


f° 1 dlnj ^ 

/ - 

J_oo Xi-r] dri 



( 100 ) 




Figure 2: (a): Normalized displacement jump associate to symmetric mechanical loading localized in iri = a = —1. Blue 
dots are COMSOL multiphysics results; (b): Normalized electric potential jump associate to symmetric mechanical loading 
localized in xi = a = 1. Blue dots are COMSOL multiphysics results. 




Figure 3: (a): Normalized shear traction associate to symmetric mechanical loading localized in = a = —1. Blue dots 

are COMSOL multiphysics results; (b): Variation of the normalized stress intensity factor with the distance a between 
te point of application of the load and the crack tip. 
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It is now possible to use (IlOQI) to obtain expressions for the stress intensity factor, and the electric 
intensity factor, 



= lim ■\/2ttxi 

xi —>-0 



( 101 ) 


Figures [2] and [31/ (&) show a comparison between the derived results and the equivalent results using 
finite element computations in COMSOL multiphysics. The point loads are assumed to be localized at 
a distance a = — 1 from the crack tip, and the normalized crack opening, electric potential jump and 
shear traction profiles are reported as functions of the spatial coordinate xi. The materials used above 
and below the crack were Barium Titanate and PZT respectively. The material parameters are quoted 
in Tabled! with those for Barium titanate obtained from Geis et al. 0 and those for PZT taken from 
Liu and Hsia 311. 


Material 

(744 (GPa) 

ei5(CM2) 

wii(C^/Nm^ 

Barium Titanate 

44 

11.4 

9.87 X 10-® 

PZT 

24.5 

14.0 

1.51 X 10-s 


Table 1: Material properties 


Good agreement between the analytical solution and the results provided by finite element analysis 
is detected in the figures for normalized crack opening, electric potential jump and shear traction field. 
The variation of the normalized stress intensity factor with the distance a bewteen the crack tip and 
the point where the loads are applied is reported in Figure/(6). 

5.2.2 Asymmetric mechanical loading 

The second example considered has point loadings at a distance a acting on the upper and lower crack 
faces. However, for this asymmetric example it is said that they both act in the same direction (see for 
example Morini et al. [l^). Mathematically this is presented as: 

(p)(xi) = 0, Ipl(xi) = . (102) 

It is important to note that, as just mentioned, most of the formulations previously proposed for fracture 
mechanics in piezoelectric bimaterials (see for example Suo et al. 0, Pak [ 3 ^, Pan Ma and Ghen 
[isl 1. do not take into account the presence of skew-symmetric components of the load acting on the crack 
faces. Consequently, in order to solve crack problems involving non-symmetric loading distribution such 
as the (IT 02 ]) . in Section [3] we have derived explicit expessions for the skew-symmetric weight functions 
matrices. 

Taking into account the loading distribution (I102L the singular integral equation (1751) becomes 

= (103) 

Using the same method as was used for the symmetric loading previously considered it is shown that 

{NA = ^ f , for - a < XI < 0 , (104) 

V M ^ Vo 1^7^34 j 

\M ) ^ V Xiy-fH34) 
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Making use of equation CH) it is possible to obtain the expression for the extended traction vector, 
t, along the interface: 


t(+)(a:i) = - 


H 


-1 /-O 


_ 1 ^[ul 

TT 7-00 xi-f] dr] 


dr], 


a 1 jj_i ( 53 H 33 


F 

TT Y iCi Xi + O I 7i/34 7 ’ 

F I 1 (S3H33H^4-lHi^' 


a;i xi + o -ff337744 — y (7 ~ ^ 3)11331134 


(106) 


It is now possible to use (110611 to obtain expressions for the stress intensity factor, K 3 , and the electric 
intensity factor, K 4 : 


= lim ^2^ f = F\[^ ^ p3i^33^44 iHlA 

\^K 4 j a:i->0 y_D j V 'Xd H 33 H 44 — y (7 ~ ^ 3 )H 33 H 34 J 


(107) 


5.2.3 Symmetric electrical loading 

We consider a symmetric system of electrical point loads on the crack faces at a distance a behind the 
crack tip. Mathematically the Dirac delta distribution is once again to represent the forces: 


Making use of equation dZSl) and the method previously used for mechanical loading gives 

5|uY“) G j a~ 1 ( 7 / 34 ^ 


dx\ 


TT Y X 1 X 1 + a \F [44 j 


When integrated, this gives 




= — arctanh — f V for — a < xi < 0 , 
TT \j a \ 7/44 / 


TT 


= —^ arctanh < / — — 


for Xi < —a. 


Xi \F[44 j ’ 

The resulting expressions for the extended traction vector is therefore 

H-i r -J—^M!-Ldr]-^ ^ ^ 


t(+Yxi) = - 


TT 7-00 xi-r] dr] 


TT y Xl Xl + o \ 1 


The stress intensity factors are then given by 


'7^3 

,7^4 


= Jim ^2Srt(+>fe)=Gy'A(“)^ 


(108) 

(109) 

( 110 ) 

( 111 ) 

( 112 ) 

(113) 


5.2.4 Asymmetric electrical loading 

Here we consider electrical loading acting in the same direction on the top and bottom crack faces at a 
distance a behind the crack tip. This loading can be written as 


(p)(xi)- 0 , |p]|(a:i) = 


0 

—2Gd{xi + a) 


(114) 
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Similarly to the case of asymmetric mechanical load (I102I1 . in order to derive the crack opening and the 
traction fields ahead of the tip correspondiiig to the awmmetric electrical loading distribution (11141) , the 
approaches proposed by Suo et al. Q , Pak , Pan Q , Ma and Chen [1^ need to be generalized to the 
case of non-symmetric electrical loading applied to the crack faces. As a consequence, the skew-symmetric 
weight function matrices derived in Section [3] are required. Assuming that the loading function is given 
by expression dm, the resulting fields obtained using equations dlSI) and (173)1 are: 


f Asf 


i9|ul|^ _ _G I a 1 f 7i/34\ 

dxi TT y xixi+ a \64H44 I ’ 

■- — arctanh ,\ , for 
TT \ a \ S4H44) 


— a < xi < 0, 


— arctanh , \ , for 

TT \ Xi yS4H44J 

t(+)fa; ) = ^ i_ ( (3-^4)i734i744 \ ^ 

^ n\xi{xi+ a){H33H44 - HI4) \S4H33H44 - 'yH^4j ' 


xi < —a. 

( {-1 - 54)H34H44 
\ 54 H 33 H 44 - ■ 

The stress intensity factor and electric intensity factor obtained are therefore 


(115) 

(116) 

(117) 

(118) 



lim i/27ra;it^“'"^(a:i) = G 

xi—>-0 


ira H33H44 - H^4 


/ (7 - 64 )H 34 H 44 \ 

\64H33H44 - lHl4j ■ 


(119) 


6 Conclusions 


A general approach for the derivation of the symmetric and skew-symmetric weight functions for plane 
interfacial cracks in anisotropic piezoelectric bimaterials have been developed. The method proposed by 
Morini et al. 171 and Pryce et al. [l^ for anisotropic elastic bodies, based on the Stroh formalism and 


Riemann-Hilbert formulation, has been extended for studying crack problems at the interface between 
dissimilar piezoelectric media. Applying this approach, explicit weight function matrices are obtained for 
an interfacial crack between two transversely isotropic piezoelectric materials, considering both the case 
where the poling direction of the two materials is perpendicular and coincident to the crack front. Since 
many poled ceramics that are commonly used in industrial applications possess transversely isotropic 
symmetry, this class of piezoelectric materials has a practical significance, and the derived weight functions 
can be used for computing the stress intensity factors corresponding to any arbitrary non-symmetric 
mechanical and electrostatic load acting on the crack faces Q . 

Symmetric and skew-symmetric weight function matrices are used together with Betti’s reciprocity 
theorem to derive integral identities relating the applied loading functions to the corresponding crack 
opening and tractions ahead of the tip. By means of the proposed method, an explicit singular integral 
formulation for the crack problem is obtained avoiding the use of Green’s function and the challenging 
numerical calculations related to such an approach. Integral identities have been derived for interfacial 
crack problems between dissimilar transversely isotropic piezoceramics subject to the two-dimensional 
state of plane strain and short circuit, having poling direction perpendicular or coincident to the crack 
front. An example of the application of the integral identities to crack problems where the loading is 
given by punctual forces acting on the faces has been performed. Explicit expressions for crack opening 
and tractions ahead of the tip corresponding to both symmetric and skew-symmetric mechanical and 
electrical loading have been obtained. The stress intensity factors associated with the introduced loading 
configurations are also evaluated. Using the derived skew-symmetric weight functions matrices, in the il¬ 
lustrative examples the effects of asymmetric loading configurations, which cannot be considered applying 
the alternative formulations available in literature, are studied. For the case where the poling direction is 
perpendicular to the 0 : 3 —axis and symmetric loading is applied at the faces, the results obtained by the 
solution of the singular integral equations are compared with those performed by finite element analysis 
using COMSOL Multiphysics. Good agreement between the analytical expression for the crack opening 
obtained by the singular integral identities and the numerical results is detected. 

Thanks to the versatility of the Stroh formalism and the generality of Betti’s theorem, the developed 
approach can be easily adapted for studying several fracture problems in piezoelectric materials with many 
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different properties without any assumption concerning the symmetry of the applied loads. Furthermore, 
the derived integral identities also have their own value from the mathematical point of view, as, to 
the authors best knowledge, such identities written in a similar explicit form for interfacial cracks in 
anisotropic piezoelectric bimaterials seem to be unknown in the literature. 
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A Derivation of the Betti formula for piezoelectric materials 


In this Appendix the derivation of the generalized Betti identity for piezoelectric materials is reported. 
Two sets of stresses, strains, electric fields and electrical displacements acting on the same physical space 
are assumed. The two sets of fields are denoted by the superscripts and respectively. The equations 
relating these fields are 24| : 

~ ^3 ^3 - ^^3 ~ ^3 ^3 ■ 

Taking the integral of (112011 over a volume, V, in combination with equation ([3]) yields 




tv 


dU = 


I ( 2 ) ( 1 )', 


{DfU^^^),3 


IV ■- 


dU. 


( 120 ) 

( 121 ) 


Making use of the Divergence Theorem and then rearranging gives 


IS 




dS' = 0, 


( 122 ) 


where S is the boundary of the volume V. 

Taking V to be a hemisphere in the upper half-plane with flat edge along the 0 : 2 —plane in equation 
(11221) leads to the following equation: 


4'>«f > + 


dxi = 0, 


'2 y ^2i “i ^2 

I X2.=0+ 

which can written in terms of the extended displacement and traction vectors used in this paper 


(123) 


/ 

Xn 


t(i).u(2)_t(2).u« dxi. (124) 

/X2=0+ ^ 

To obtain equation (E51) and are taken as the physical fields and then = RU and = RS. 
Equation dia is derived using a semi-circular domain in the lower half-plane. 
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B Transversely isotropic piezoelectric materials: explicit ex¬ 
pressions for matrices Q, R and T 

In this Appendix, explicit expressions for the stiffness, permittivity and piezoelectric tensors correspond¬ 
ing to transversely isotropic piezoelectric materials with poling direction parallel to the X 2 and axes are 
reported. Using these tensors, general forms for the Stroh’s eigenvalues matrices, the surface admittance 
tensor and the bimaterial matrices H and W are obtained. 


B.l Poling direction parallel to the X 2 —axis 

When considering transverse isotropic materials with poling direction parallel to the a; 2 -axis the stiffness 
tensor, C, simplihes to 


/Cii 

C'12 

Ci3 

0 

0 

0 \ 

C'12 

C22 

C'12 

0 

0 

0 

C'ls 

C'12 

Cii 

0 

0 

0 

0 

0 

0 

C44 

0 

0 

0 

0 

0 

0 

(Cii - ^13)72 

0 

V 0 

0 

0 

0 

0 

C'44/ 


and the permittivity and piezoelectric tensors simplify to 



^UJii 

0 


(0 

0 

0 

0 

0 

ei 6 \ 

UJ = 

0 

W22 

0 , e = 

621 

622 

621 

0 

0 



1 0 

0 

UJllJ 


0 

0 

616 

0 

0/ 


This is the same system as used in Hwu 
reduce to 


22 |. 


Using these conditions the matrices from equation m 


fCn 

0 

0 

0 ^ 


/ 0 

C '12 

0 

62i\ 

0 

C '44 

0 

616 

, R = 

C '44 

0 

0 

0 

0 

0 

(Cii - Ci3)/2 

0 


0 

0 

0 

0 

V 0 

616 

0 

-Wil/ 


\ei6 

0 

0 

0 / 


/c'44 00 0 \ 

0 C22 0 622 

0 0 C '44 0 

\ 0 622 0 —UJ 22 ) 


Note that in this case the poling direction is perpendicular to the crack plane reported in Fig. [TJ The 
antiplane component can clearly be decoupled from the rest of the elasticity components and all effects 
caused by the electric field. This means that the Mode III tractions and displacement have no dependency 
on the electric field and therefore behave in the same way as they would in an elastic material with no 
piezoelectric effect. 

Only the in-plane components and electrical effects are considered. That is: u = {ui,U 2 ,(j>)'^ £ind 
t = {<J 2 i,cr 22 , 02 )^■ The decoupled part of the eigenvalue problem ([5]) now has matrices 


Q 


/Cii 

0 

0 ^ 



( ^ 

C'12 

0 

C44 

616 

, R 

= 

Cm 

0 

V 0 

616 

-Wll/ 



\e2i 

0 



(Cm 

0 

0 




T = 

0 

Cm 

622 




1 0 

622 

—W22/ 




(125) 


The extended Stroh formalism described in Section 2 of the paper is an effective tool for finding 
an expression for the surface admittance tensor, B. However, the eigenvalue problem obtained is not 
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always straightforward to solve. For non-piezoelectric materials an alternative approach is the Lekhnitskii 
formalism introduced by Lekhnitskii which gives a specific normalisation of the eigenvalues obtained 
from the eigenvalue problem (jS]) 33| . The Lekhnitskii formalism was extended to the piezoelectric setting 
by Hwu 23 where it was used to find the surface admittance tensor for the poling direction parallel to 
the a; 2 -axis 



/ 

Bn 

iBi 2 

iBi4 

B = 


-iBi 2 

B22 

B24 


V 

-iBi4 

B24 

B 44 


(126) 


The expressions for the components of the matrix B found by Hwu are quoted in the Appendix C. 

With an expression for B it is now possible to construct the bimaterial matrices required. The 
bimaterial matrix H can be written as 



( Hn 

iHi2 

iHi4 

H = 

-iHi2 

H 22 

H 24 


\^-iHi4 

H 24 

H 44 

H(yoL - [-^act]/ [^Cta]//7 

for a 

Hip = [Bi 

p]i - [Bi 


for /3 


where 


H 2 A = [H 24 ]/ + [ 524 ]//. 

The matrix W has the same structure as H, that is 


W = 


/ Wn iWi2 iWii'^ 
—iWi 2 W22 IF24 

y— iWi4 W24 W44 


where 


kFcKQ, — for cr — 1,2,4, 

Wip = [B 4 p]i + [B 4 p]ii, for /3 = 2,4, 
W 24 = [S 24 ]/ — [.B24]//- 


(127) 


(128) 


B.2 Poling direction parallel to the X 3 —axis 

When considering transverse isotropic materials with poling direction parallel to the xs-axis the stiffness 
tensor, C, simplifies to 


(Cn 

C12 

Ci3 

0 

0 

0 \ 

C 12 

Cn 

Ci3 

0 

0 

0 

C'la 

Ci3 

C33 

0 

0 

0 

0 

0 

0 

C44 

0 

0 

0 

0 

0 

0 

C44 

0 

V 0 

0 

0 

0 

0 

(Cn - Ci2)/2y 


and the permittivity and piezoelectric tensors simplify to 



^un 

0 


/o 

0 

0 

0 

ei5 

UJ = 

0 

Wn 

0 , e = 

0 

0 

0 

ei5 

0 


1 0 

0 

W33/ 

\e3i 

631 

633 

0 

0 
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Using these conditions the matrices from equation d®]) reduce to 


/Cn 

0 

0 

^ ) 


/ 

0 

Ci 2 

0 

o\ 

0 

(Cii - C12 

)/2 0 

0 

, R = 

(Cii- 

- Ci2)/2 

0 

0 

0 

0 

0 

C44 

ei 5 



0 

0 

0 

0 

V 0 

0 

ei 5 

-UJii/ 


l 

0 

0 

0 

0/ 



^{Cii — 

Ci2)/2 

0 0 

0 

\ 





^ ^ 0 Cii 0 0 

0 0 C44 ei5 

\ 0 0 ei5 —oJiiJ 

Note that in this case the poling direction coincides with the front of the crack reported in Fig. [TJ For 
this particular case it is possible to decouple the Mode I and Mode II components of the displacement and 
stress helds from the Mode III fields and piezoelectric effects on the material. This means that the Mode 
I and II helds will behave similarly to those for purely elastic materials with no piezoelectric behaviour. 

Only the out-of-plane and piezoelectric components are considered. In this case Q, R and T are 
reduced to 2 x 2 matrices: 

Q = T = , R = 0 . (129) 

ycis —ujiij 

As a consequence, for this case the surface admittance tensor, B assumes the form 


/ B33 B34 

yR34 B44 


Explicit expressions for the components of B are given in the Appendix C. 
The bimaterial matrices H and W can be computed and have the form 


(H33 H34\ ^^(^3H33 1H34\ 

\H34 H44)' \lH34 S4H44)' 


The components of these matrices are given by 


Ba^ — 

6 a = 


[Bal3]l + [Bal3]ll, for o, /3 = 3,4, 
[Baa]l — [Bo 


II 


[Baa]l + [Baa]ll ’ 


for a = 3,4, 


7 = 


[. 634 ]/ — [B34]ii 
H 3i 


(130) 


(131) 


C Explicit expressions for B = zAL ^ 

In this section we give the expressions for the surface admittance tensor B corresponding to the two- 
dimensional state of plane strain and short circuit (£3 = 0, E 3 = 0), as seen in Hwu [l^. We only give 
the expressions for the decoupled part of the tensor which contains the piezoelectric behaviour of the 
material. 


C.l Poling direction parallel to the a: 2 -axis 

As stated in Section 2, the general form of the matrix B = for transverse isotropic materials with 

poling direction parallel to the X 2 -axis is 



/ 

Bii 

iBi 2 

iBi 4 

B = 


-iBi 2 

B22 

B24 


V 

-iBi 4 

B24 

B44 


( 132 ) 
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Here, explicit expressions for the components of this matrix (found in Hwu are given. 

The following components of the compliance tensor: S, piezoelectric strain/voltage tensor: g, and 
dielectric non-permittivity: /3, are introduced: 


C22 (e2iC'ii — 6221^12)^ 

C11C22 - Cl^ ~ C*[CiiC22 - cl^] ’ 
621622 [C*!! — C'12] W22C'i2 

~ C*[CiiC 22 -Cl^] c* ’ 

621 [^11 ~ ^12] I <^22Cll 

C*[CnC22 - 0 ^ 2 ] C* ’ 


S'66 = 

9 '21 — 
9 '22 = 


0^11 

6i6 + ^11^44 

621(^11 ~ 6221^12 

c* 

6221^11 — 62 iC'i 2 


5^16 ~ 
/3'ii = 
P'22 = 


616 

^16 + W 11 C 44 

C 44 

®16 ‘^11^44 ’ 

gllC22 - Cl^2 
C* 


where 

C* = (621 + 622 )C'i 1 — 2e2l622C'l2 + W 22 [C'llC '22 — <^ 12 ]- 

Through using the Lekhnitskii formalism extended to piezoelectric materials the eigenvalues, p, are 
found through the equation 

kp 2 - m 3 = 0, (133) 

where k, p 2 and m 3 are functions of p and are given by 

k = S'iip'^ + {2S '12 +SkojP^ + S' 22 , rns = -{ 9 ' 21 + 9 'w)P^- 9 ' 22 ^ P 2 =-{P'nP^ + P' 22 )- (134) 


This sextic equation must be solved numerically but is easily shown to have roots of the form 


P2=a2+ iP2, P3 = -02 + iP2, Pi = iPi- 


(135) 


With the eigenvalues known, Hwu 
B. It was shown that 


22 | proceeded to find explicit expressions for the components 


Bii = 2S'iil-m{plij2 + {plri2 - plm)}/\ 

B 22 = 2Im{[72P2P4?74 + ( 72 P 4 - liP2)P2\/P2Pi}/\ 

B 44 = -2P'^^lui{p2P2'r]2 + P 2 Pi{m - Pi)}/\ 

B 24 = 2 / 3 ' 33 lm{p 2 P 2 » 72»74 + P 2 Pi'lf 2 {r ]2 - V4)}/\ 

Bi 2 = S'i 2 + 2 Re{[ 72 P 4??2 + ( 74 P 272 - l 2 Pi'ni)\/P 2 Pi}/\ 
Bli = -g'le + 2P'j^^Re{p 2 P 2 V 2 - V2V4iP2 - P4)}/\ 

where 

A = 2Re{p2?72 + {P4P2 - P 274 )}, 

Ik = S '22 + 9 ' 22 Vk, Vk = for k=2,4. 

m3[pk) 


25 
















C.2 Poling direction parallel to the xs-axis 

The decoupled part of the surface admittance tensor, B = containing Mode III fields and the 

piezoelectric effect for transverse isotropic materials with poling direction parallel to the xa-axis is: 


/ Baa S34 

\yB34 B44 


Explicit expressions for the components of this matrix are given by: 


B33 = 


UJll 

®15 + <^ 111^44 ’ 


B44 

B34 


-C44 

®15 + <^ 111^44 
ei5 

®15 + <^ 111^44 


(136) 


D Inverse Fourier transforms 


In order to derive explicit expressions for equations dlol) and m, we need to compute the inverse Fourier 
trasform of terms of the form *sign(^)/(^), |CI/(C) and Using the Fourier convolution theorem 

and the knowledge that the inverse Fourier transform of sign(^) is given by —i/{'KXi) gives: 

B"^[isign(0/(^)] = iJ-^[sign(^)] 

= i 

1 


-) 

TTXi / 




r fiv) 

Loo Xi-rj 


dry. 


(137) 


The inverse Fourier transform of |C|/(0 is found: 


-^-MlCI/( 0 ] = -^-Msign(C)] *-F-MmC)], 

=(—)#■ 

axi 

Finally, the inverse Fourier transform of zsign(^)|^|/(^) is given by 

J^"^[rsign(^)|^|/(C)] = B“^[rC/(C)], 

= -LL 

dxi' 


(138) 


(139) 


E Expressions for matrices M and N 

In this Appendix explicit expressions for the matrices M and N are quoted. They have the form 

M = ^ (M' + rsign(OM"), (140) 

N = M (N' + rsign(ON"), (141) 

where 

D = HI^H22 + Hl^Hn + - HiiH22H44 - 2 H 14 H 42 H 24 . (142) 


26 






The matrices M',M",N' and N" have the form 


where 



/Mil 

0 

0 \ 

( 

Mi2 

Mi4\ 

M' = 

0 

M22 

II 

M 21 

0 



1 0 

M42 

M44) 

\M 41 

0 

0 / 



/iVii 

0 

0 \ 

( ^ 

-Y 12 

-Yi4\ 

N' = 

0 

N 22 

II 

-iVi2 

0 



1 0 

-Y 24 

7 V 44 / 

A^14 

0 

1 / 


Mil = Vhii(M |4 - H 22 H^i) + Wi 2 {Hi 2 Hii - M14M24) - Wu{Hi 2 H 2 A - H14H22), 
M22 = Wi2{Hi2Hii - M14M24) - Vh 22 (i?lli ?44 - Hli) - W 2 i{Hi 4 Hi 2 " H11H24), 
M44 = Wu{HuH22 - H12H24) - W 24 {HuHi 2 - i^ll^^ 24 ) " Wii{H^iH22 " HI2), 
M24 = Wu{HuH 24 - H12HA4) + W 24 {HiiHii - + W 44 {HuHi 2 - H11H24), 

M42 = Wi 2 {Hi 2 H 24 ~ H14H22) + W 22 {HhHx 2 ~ H11H24) + ^^24(^11^22 — H12), 


(143) 


(144) 


(145) 

(146) 

(147) 

(148) 

(149) 


Mi 2 = Wi 2{H22H44 - H^i) - W22{Hi2H44 - H 14 H 24 ) + W24{Hi2H24 - H 14 H 22 ), (150) 

Mi 4 = WiiiHli - H 22 H 44 ) + W24{Hi2H44 - H 14 H 24 ) - W 44 (i 7 l 2 i 724 - HUH 22 ), (151) 

M 21 = Wn{Hi 2 H 44 - H 14 H 24 ) - Wi 2 {HiiH 44 - HI 4 ) - Wi 4 {HuHi 2 - HnH24), (152) 

M 41 = Wn{Hi2H24 - H 14 H 22 ) + Wi2{Hi4Hi2 - H 11 H 24 ) + Wu{HiiH 22 - HI 2 ), (153) 


IVii = H 22 H 44 - HI 4 , (154) 

N 22 = H 41 H 44 - Hl^, (155) 

IV 44 = i7iii722 - HI 2 , (156) 

N24 = H4lH24-H44Hi2, (157) 


lVi 2 — H 12 H 44 — H 14 H 24 , 
Ni4 = H 12 H 24 — Hi 4 H 22 - 


F Expressions for matrices Y and Z(^) 

Explicit expressions for the matrices Y and Z(^) are given by: 

Y = -U~^W = 1 (hH33H44 - IHI 4 77447734(7 - ^4) \ 

2 2(77337744 — 7714 ) y 77337734(7 — ^ 3 ) 54 H 33 H 44 — ^H^ 4 j 


m 


lel ( H 44 -H 34 \ 

H 33 H 44 — 77I4 7734 H 33 J 


(158) 

(159) 


(160) 

(161) 
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